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The relativistic coupled-cluster single-double method is used to calculate the dependence of fre- 
quencies of strong i?l-transitions in many monovalent atoms and ions on the fine-structure constant 
a. These transitions are used in the search for manifestations of the variation of the fine-structure 
constant in quasar absorption spectra. Results of the present calculations are in good agreement 
with previous calculations but are more accurate. 
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I. INTRODUCTION 

Theories unifying gravity with other interactions sug- 
gest a possibility of temporal and spatial variations of the 
fundamental constants of nature; a review of these theo- 
ries and results of measurement can be found in Ref. [l[ . 
A very sensitive many-multiplet (MM) method to search 



for the variation of the fine-structure constant a 



by comparison of quasar absorption spectra with labora- 
tory spectra has been suggested in Refs. 0, 0]. 

Using this method, strong evidence that the fine- 
structure constant might be smaller about ten billion 
years ago was found [3, 0, S H, 0, [!]■ This result was 
obtained from an analysis of data from the Keck tele- 
scope in Hawaii by the group based at the University of 
New South Wales in Australia. However, an analysis of 
data from the VLT telescope in Chile, performed by dif- 
ferent groups [9, l'^ using the same MM method, gave a 
null result. There is an outgoing debate in the literature 
about possible reasons for the disagreement. 

The MM method requires calculation of relativistic 
corrections to frequencies of atomic transitions to reveal 
their dependence on the fine-structure constant. All cal- 
culations used in the analysis so far were performed by 
a single research grorip based at the University of New 
South Wales [3, M, [M 113, d, E, 13 ■ Owing to the 
importance of detecting variations of fundamental con- 
stants and the disagreement mentioned above, it is im- 
portant that the related atomic calculations be verified 
independently. 

A positive development in this direction is a recent in- 
dependent calculation of the relativistic energy shifts in 
the ion Fe II [2l[ . Fe II is the single most important ele- 
ment in the analysis of quasar absorption spectra. It has 
lines which move in opposite directions when a varies; 
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moreover, the shifts in Fe II are relatively large. In prin- 
ciple, subject to sufficient statistics, Fe II alone could 
serve as a probe of variation of the fine-structure constant 
in quasar absorption spectra [2l| . However, calculations 
for Fe II are difficult due to the large number of valence 
electrons. 



In a recent work [23| independent calculations of the 
relativistic energy shifts for monovalent and diavalent 
atoms of astrophysical interest were reported. This work 
also presented a detailed analysis of Breit contributions 
to the relativistic energy shift. Note that, despite some 
overlap of authors of early [11 [H, [H, [II M, El, HI 
and recent [U, [S^l works, the later can be regarded as 
independent since they use completely independent sets 
of computer codes and different methods of calculation. 



In present work, we study the role of higher-order cor- 
relations in relativistic energy shifts of monovalent atoms 
and ions of astrophysical interest. We use a linearized 
coupled-cluster method in the single-double approxima- 
tion and third-order many-body perturbation theory to 
perform the calculations. We demonstrate that including 
higher-order correlations significantly improves the accu- 
racy of the transition energies. The values of relativistic 
energy shifts found in the present work are in good agree- 
ment with previous calculations but are more accurate. 
Note also that the present work can be considered as an 
independent verification of earlier calculations. This is 
because the calculations are performed using a method 
that has never before been used for this purpose. 



Apart from the important task of calculating relativis- 
tic energy shifts of atomic frequencies of astrophysical 
interest, this work can be considered as another demon- 
stration of the power of the single-double method. It 
applies the method for the first time to the ions C IV 
and Zn II. 
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II. METHOD 

It is convenient to present the dependence of atomic 
frequencies on the fine-structure constant a in the vicin- 
ity of its physical value ag in the form 



iLj{x) = luq + qx, 



(1) 



where wq is the laboratory value of the frequency and 
X = (a/ao)^ ~ Ij 9 is the coefficient which is to be found 
from atomic calculations. Note that 



dhJ 

dx 



x=0 



To calculate this derivative numerically we use 



25 



(2) 



Here 6 must be small to exclude non-linear in terms. 
In the present calculations we use 6 = 0.05, which leads 
to 



q^W (w(+0.05) - w(-0.05)) . 



(3) 



To calculate the coefficients q using ^ , a must be varied 
in our computer codes. Therefore, it is convenient to use 
a form of the single electron wave function in which the 
dependence on a is explicitly shown (we use atomic units 
in which e = h = l,a ^ 1/c) 



njlm 



(4) 



where n is the principal quantum number and an index 
V replaces the three- number set n,j,l. This leads to a 
form of radial equation for single-electron orbitals which 
also explicitly depends on a: 



dfv 



dr r 

dgy K 

dr 



+ —fv{r)^ 2 + a'{€,-VHF) 9v{r) = 



+ -V>/F)/.(r) = 0, 



(5) 



here k = + 1/2), and Vhf is the H artree- 

Fock potential. Equation (O with a = a^^JS -\- 1 is used 
to construct a full set of single-electron orbitals. 

As a first step, equation © is used to calculate self- 
consistently single-electron states of the closed-shell core. 
Then, this equation is used to calculate a complete set 
of B-spline single-electron basis orbitals. We use 40 B- 
splines of order fc = 9 in a cavity of radius 4Ga_B with 
angular momentum I up to Imax = 5. More details on 
the use of B-splines in atomic calculations can be found 
in Ref. 

This basis set is used to perform calculations with 
the linearized couple-cluster single-double method sup- 
plemented by third-order many-body perturbation the- 
ory. The all-order single-double (SD) method was dis- 
cussed before in Refs. [H [H, [ll, 113, HI] and third- 
order many-body perturbation theory was discussed in 



Refs. [1^, . We also include results from second-order 
MBPT to uncover the role of higher-order correlations 
and to estimate the numerical uncertainty due to corre- 
lations. 

The SD equations are written for the coefficients of the 
expansion of the many-electron wave function in terms of 
single and double excitations from the reference Hartree- 
Fock wave function. They contain Coulomb integrals be- 
tween single-electron basis orbitals as parameters. The 
equations are first solved iteratively for the atomic core. 
When convergence is achieved, a similar procedure for 
valence states of interest is performed. 

Many-body perturbation theory is used to calculate 
the third-order diagrams missed in the SD method 
(E3extra). Therefore, the present calculations are com- 
plete through third order while selected classes of higher- 
order diagrams are included in all orders. A number of 
earlier calculations (see, e.g. [H [H, [H, [13, [H, [11] ) prove 
that this approach gives very good accuracy for a wide 
range of monovalent atoms and ions. 

The SD equations depend on the fine-structure con- 
stant only implicitly, via the values of the Coulomb in- 
tegrals, which are affected by the change of the single- 
electron basis states @ due to change of a in the Hamil- 
tonian (O. In other words, changes in a lead to changes 
in the basis orbitals, while the SD equations remain un- 
changed. The same is true for terms in the MBPT ex- 
pansion. However, the change of basis means that all 
SD and E3 calculations must be redone. The E3 calcu- 
lations must be repeated from scratch, but the SD iter- 
ations can be restarted from previous calculations for a 
different value of a. Since the change of a is small, [see 
Eq. ^] only a few iterations are needed to solve the SD 
equations for new value of a. This significantly speeds 
up the calculations. 

Although the method used in present work has been 
used many times before, the set of computer codes used 
in present work is new and independent of previous ver- 
sions. The main reason for developing new codes was the 
need to have them in a form that allows easy modifica- 
tions, as e.g. the change of the fine-structure constant 
in this calculations, or inserting extra operators (Breit 
interaction, specific mass shift operator, etc.) for future 
projects. Some features of the present realization of the 
methods are presented in the appendix. For third-order 
MBPT calculations, we apply a universal algorithm that 
uses a numerical description of the E3 MBPT diagrams 
and the same piece of computer code to calculate all of 
them. This algorithm will be described in detail else- 
where I30i. 



III. RESULTS 

Results of calculations of energy levels of C IV, Na I, 
Mg II, Al HI, Si IV, Ca II and Zn II are presented 
in Table [D Removal energies of the lowest s and p- 
states, which are important for the analysis of quasar 
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TABLE I: Removal energies of the lowest s and p states of C IV, Na I, Mg II, Al III, Si IV, Ca II and Zn II in different 
approximations, comparison with experiment (cm^^). 



Atom 


Z 


State 


RHF 


E2 


A" 


SD 


A" 


SD+E3 


A" 


Expt* 


C IV 


6 


2s 


519255 


520189 


11 


520233 


55 


520231 


53 


520178 






2pi/2 


454054 


455640 


-54 


455732 


38 


455727 


33 


455694 






2p3/2 


453927 


455509 


-78 


455600 


13 


455595 


8 


455587 


Nal 


11 


3s 


39952 


41229 


-220 


41437 


-12 


41450 


1 


41449 






3pi/2 


24030 


24412 


-81 


24486 


-7 


24489 


-4 


24493 






3P3/2 


24014 


24394 


-82 


24468 


-8 


24471 


-5 


24476 


Mg II 


12 


3s 


118824 


121076 


-192 


121278 


10 


121273 


5 


121268 






3pi/2 


84294 


85453 


-145 


85585 


-13 


85586 


-12 


85598 






3p3/2 


84204 


85357 


-150 


85488 


-19 


85489 


-18 


85507 


Al III 


13 


3s 


226396 


229319 


-127 


229489 


43 


229464 


18 


229446 






3pi/2 


173687 


175608 


-155 


175751 


-12 


175741 


-22 


175763 






3p3/2 


173452 


175361 


-168 


175504 


-25 


175494 


-35 


175529 


Si IV 


14 


3s 


360614 


364033 


-60 


364172 


79 


364132 


39 


364093 






3pi/2 






1 

-LOU 




-4 




9fi 








3p3/2 


289606 


292183 


-161 


292317 


-27 


292296 


-48 


292344 


Ca II 


20 


4s 


91440 


96173 


425 


96097 


349 


95577 


-171 


95748 






4pi/2 


68037 


70680 


123 


70761 


204 


70491 


-66 


70557 






4p3/2 


67837 


70449 


115 


70529 


195 


70262 


-72 


70334 


Zn II 


30 


4s 


135134 


143835 


-1055 


144618 


-272 


145334 


444 


144890 






4pi/2 


90524 


95249 


-1161 


96184 


-226 


96613 


203 


96410 






4p3/2 


89787 


94372 


-1164 


95311 


-225 


95728 


192 


95536 



"A — -Bcalc — -Eexpt 

''NIST, Ref. [H 



absorption data, are presented in different approxima- 
tions. These include the relativistic Hartree-Fock ap- 
proximation (RHF), second-order many body perturba- 
tion theory (E2), single-double (SD) approximation, and 
SD supplemented by third-order many-body perturba- 
tion theory (SD-I-E3). For each approximation the dif- 
ference between theoretical and experimental energies is 
presented in the columns headed A. 

The Mg II, Al III, and Si IV ions represent an isoelec- 
tronic sequence of sodium. This sequence was considered 
in detail in Ref. [1^ . The results of present work are in 
good agreement with previous calculations. Some small 
difference can be attributed to the difference in numer- 
ical procedures and numerical parameters (such as the 
number of splines, the cavity radius, maximum angular 
momentum, etc.) 

The results in Table |I] show that correlations are large 
and are strongly dominated by second-order MBPT. 
However, inclusion of higher-order correlations is impor- 
tant and leads to a significant reduction in the differences 
between theoretical and experimental energies. 

Table |lT] presents relativistic energy shifts {q coeffi- 
cients) in the same approximations as the energies in 
Table H] One can see that the role of correlations is 
much less important for the q-coefRcients than for en- 
ergies. While second-order correlations give some rather 
small contribution to q, contributions from higher-order 
correlations are practically negligible in most cases. Only 
for Zn II the higher-order contributions are significantly 
larger than the uncertainty of the calculations. Thus, 



TABLE II: Relativistic energy shifts (g-coefBcients, see 
Eq. (m for the lowest s and p states of C IV, Na I, Mg II, 
Al III, Si IV, Ca II and Zn II in different approximations 
(cm-i). 



Atom 


Z 


State 


RHF 


E2 


SD 


SD-I-E3 


C IV 


6 


2s 


243 


244 


244 


244 






2pi/2 


139 


141 


142 


142 






2p3/2 


11 


10 


10 


10 


Na I 


11 


3s 


51 


57 


59 


59 






3pi/2 


12 


13 


13 


13 






3p3/2 


-5 


-5 


-6 


-5 


Mg II 


12 


3s 


181 


193 


194 


194 






3pi/2 


70 


74 


74 


74 






3P3/2 


-20 


-23 


-23 


-23 


Al III 


13 


3s 


405 


421 


422 


422 






3pi/2 


197 


203 


203 


203 






3p3/2 


-39 


-45 


-45 


-44 


Si IV 


14 


3s 


753 


773 


774 


774 






3pi/2 


416 


425 


425 


425 






3p3/2 


-54 


-62 


-62 


-61 


Ca II 


20 


4s 


354 


396 


394 


392 






4pi/2 


161 


176 


176 


175 






4p3/2 


-41 


-58 


-59 


-56 


Zn II 


30 


4s 


2352 


2863 


2873 


2910 






4pi/2 


990 


1271 


1333 


1359 






4p3/2 


224 


354 


412 


427 



calculations of the g-coefficients are more stable than the 
calculations of energies and, consequently, more accurate. 
Breit corrections to the frequencies of the transitions 
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TABLE III: Contributions to the relativistic energy shifts (g- 
coefficients, see Eq. ([1} for the s - p transitions in C IV, Na I, 
Mg II, Al III, Si IV, Ca II and Zn II (cm-^) 



Atom 


z 


Transition 


Dirac" 


Breit' 


Total 


C IV 


5 


2« — 2t)i /o 


102 


13 


115 








234 


-12 


222 


Na I 


11 


3,'i 3t)i /o 


46 


-1 


45 








64 


_2 


62 


Mg II 


12 


3s — 3pi/2 


120 


1 


121 






3s - 3p3/2 


217 


-5 


212 


Al III 


13 


3s - 3pi/2 


219 


5 


224 






3s " 3p3/2 


467 


-9 


458 


Si IV 


14 


3s - 3pi/2 


348 


13 


361 






3s - 3p3/2 


835 


-12 


823 


Ca II 


20 


4s - 4pi/2 


218 


3 


222 






4s - 4p3/2 


450 


-4 


446 


Zn II 


30 


4s - 4pi/2 


1546 


-5 


1541 






4s - 4p3/2 


2472 


-20 


2452 



"This work, see Table |lT] 
''Savukov and Dzuba, Ref. 



considered in present paper were calculated in Ref. [22] ■ 
These corrections are larger than the uncertainty of 
the g-coefficients from omitted higher-order correlations. 
Therefore, they must be included for accurate results. 
Table IIIII summarizes all significant contributions. The 
Dirac contributions in this Table are based on averaging 
the SD and SD-I-E3 approximations from Table [Til Breit 
corrections are taken from Ref. [22| . 

Our final results are presented in Table IIVI An esti- 
mate of the numerical uncertainty is also given and re- 
sults are compared with previous calculations. There are 
two sources of numerical uncertainty. One is omission of 
certain higher-order correlations and numerical accuracy 
of the SD and SD+E3 calculations. This uncertainty 
was estimated by comparing the results in the SD and 
SD+E3 approximations. Another source of uncertainty 
is the accuracy of calculation of Breit contribution. Breit 
contributions are calculated very accurately within the 
relativistic Hartree-Fock approximation (see Ref. [l^l for 
details). The only uncertainty which may come from 
the Breit interaction is due to the fact that correlation 
corrections to the Breit interaction are ignored. How- 
ever, these contributions are small for the relatively light 
atoms considered in present paper. A strong argument 
that the Breit contribution is calculated sufficiently ac- 
curately was presented in Ref. (2^ , where it was demon- 
strated that the inclusion of the Breit interaction brings 
the theoretical fine-structure intervals into perfect agree- 
ment with experiment. We use a rather conservative esti- 
mate of 20% for the accuracy of the calculations of Breit 
contributions. 

Table [IV] also presents results from previous calcula- 
tions. All previous calculations were done in second-order 
of MBPT. The work of Ref. [H also included Breit cor- 
rections. Results of the present paper are in excellent 
agreement with this work. Some difference for Zn II is 



TABLE IV: Relativistic energy shifts (g-coefficients, see 
Eq. m for the s - p transitions in C IV, Na I, Mg II, Al III, 
Si IV, Ca II and Zn II (cm~^); comparison with other calcu- 
lations. 



Atom 


Z 


Transition 


This 


Other 










work 


Savukov" 


Berengut'' 


C IV 


6 


2s 


— 2pi/2 


115(2) 


115 


104(20) 






2s 


— 2p3/2 


222(2) 


221 


232(20) 


Na I 


11 


3s 


- 3pi/2 


45(0) 


44 


45(4) 






3s 


- 3p3/2 


62(0) 


61 


63(4) 


Mg II 


12 


3s 


- 3pi/2 


121(1) 


120 


120(10) 






3s 


- 3p3/2 


212(1) 


211 


211(10) 


Al III 


13 


3s 


- 3pi/2 


224(1) 


223 


216(14) 






3s 


- 3p3/2 


458(2) 


457 


464(30) 


Si IV 


14 


3s 


- 3pi/2 


361(2) 


360 


346 






3s 


- 3p3/2 


823(2) 


823 


862 


Ca II 


20 


4s 


- 4pi/2 


222(1) 


222 


224 






4s 


- 4p3/2 


446(3) 


450 


452 


Zn II 


30 


4s 


- 4pi/2 


1541(7) 


1585 


1584(25) 






4s 


- 4p3/2 


2452(13) 


2488 


2479(25) 



"Second-order-l-Breit, Ref. [23| 

"Second-order, compilation of previous results presented in 
Ref. [531 



due to higher-order correlations. The results are also 
in good agreement with early calculations compiled to- 
gether in Ref. [13] . The principal sources of the small dif- 
ferences between the present results and those of Ref. [2^ 
are Breit and higher-order correlation corrections. 



IV. CONCLUSIONS 

In the present work, we studied the role of higher- 
order correlations on relativistic energy shifts of atomic 
frequencies used in the search for variations of the fine- 
structure constant in quasar absorption spectra. We have 
demonstrated that the higher-order correlations are im- 
portant for energies, bringing theoretical values into bet- 
ter agreement with experiment. However, higher-order 
correlations give very small contributions to the relativis- 
tic energy shifts (g-coefficients) in all cases except Zn II. 
Results of the present work are in good agreement with 
previous calculations but are more accurate. 
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APPENDIX A: EFFICIENT WAY OF 
CALCULATING THE SD AND MBPT TERMS 

The coupled-cluster SD method combined with MBPT 
has led to an accurate description of many properties of 
monovalent atoms; however, it is very demanding com- 
putationally. For example, for Zn II the total number of 
single-electron basis states used in present calculations is 
319. The total number of non-zero distinctive Coulomb 
integrals is about 1.7 x 10*. All of them are used in both, 
the SD equations and the MBPT expansion. Calculating 
all Coulomb integrals in advance and keeping them in 
computer memory is practically impossible due to huge 
demand for computer memory. On the other hand, cal- 
culation of Coulomb integrals from single-electron basis 
functions every time they are needed makes the calcula- 
tions unacceptably slow. This is even more so in the case 
of the relativistic energy shifts considered in present pa- 
per, since we need to run all relevant codes several times 
for several different values of the fine-structure constant 
a. 

To improve the efficiency of the codes, we use an ap- 
proach in which Hartree screening functions Y rather 
than Coulomb integrals are calculated in advance and 
kept in memory for efficient calculation of Coulomb in- 
tegrals (a similar approach was used in Ref. [l^l)- The 
Hartree screening function Y is defined as 

y^nmir) = y ^(/„(r')/™(r') 

^o?g^(r')g^{r'))dr\ (Al) 



where r< = min(r, r') and r> = max(r, r'). We also need 
the p functions: 

P,i{r) = /j(r)//(r) + (x^ g,{r)gi(r). (A2) 



Our typical coordinate grid consists of about 1000 
points. Usually all points are used to calculate Y-^nmij'^- 
However, there is no need to keep all points for succes- 
sive calculations of the Coulomb integrals. It turns out 
that very little lose of accuracy is caused by using a sub- 
set of points defined as every 4th point in the interval 
1/Z < r < i?cavity, whcrc i?cavity IS the radius of the cav- 
ity in which the i?-splines basis orbitals are defined. By 
cutting off points at short distances and using only every 
4th point in between we reduce the number of points by 
an order of magnitude. Then, the Coulomb integrals are 
calculated in an extremely efficient way as 

qkijlmn) = ^ Pji{ri)Yk,nn{ri)wi. (A3) 

i=l 

Here fj, « 100 is number of points on the sub-grid and 
Wi are weight coefficients corresponding to a particular 
method of numerical integration. Note that only one 
of two integrations for Coulomb integrals is done on a 
reduced sub-grid. The initial integration (jAip is done 
using all grid points. 
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